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Abstract 

In this communication we show that the surface tension of the real fluids of the Lennard-Jones 
type can be obtained from the surface tension of the lattice gas (Ising model) on the basis of the 
global isomorphism approach developed earlier for the bulk properties. 
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The lattice models play the important role in the understanding the behavior of the 
real systems which are too complex to be treated in controllable way. In general the direct 
comparison between real continuum systems and their discrete models is impossible. That is 
why the situations where one may to establish the quantitative relations between the physical 
properties of these systems are of great importance. Recently, the conception of the global 
isomorphism between Lennard- Jones fluids and the short-ranged Ising model (lattice gas) 
was proposed in [I]. It is based on the simple geometrical reformulation in [2j the results 
of the profound analysis of the thermodynamic data for real and model systems made by 
E. Apfelbaum and V. Vorob'ev in a series of papers [3j H]. These authors stressed the 
necessity to pay special attention on the linear regularities characteristic for the molecular 
fluids known more than a century. These linearities are the rectilinear diameter law [5] 
and, especially, the Zeno line linearity [6HH] in the simplified form known as the Batchinsky 
law [9\. In [2] it was shown how these phenomenology could be casted in clear geometrical 
picture via the mapping between the set of the states of the fluid and the lattice gas of the 
form: 

The objective of this paper is to use the global isomorphism transformation to relate the 
surface tension of the liquid- vapor interface of the lattice gas (Ising model) with that of the 
LJ fluid. The results obtained for the binodal in [10] give the ground to use the approach 
for the surface tension since from the thermodynamical point of view this physical quantity 
is nothing but the thermodynamic potential of the surface. It has obvious meaning in the 
lattice gas model and due to the existence of the exact solutions it is possible to introduce 
the thickness of the surface. 

The following mapping representing the global isomorphism between the thermodynamic 
states of the LG and the LJ fluid was constructed: 
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Here n and T are the density and the temperature of the fluid, t is the temperature variable 
of the LG normalized to the critical temperature so that t c — 1. We also use the standard 
dimensionless values for T and n of the LJ fluid [H]. T* and are the parameters of the 
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linear zeno-element: 

The coordinates of the CP for the fluid are: 

2(1 + a ) I + a 



^A/" — 1Y max "•" iV l "T • • • 

where A/* = mi x 7712 is the size of the lattice. The thermodynamic potential is augmented 
with the surface contribution is: 

VT In Sy(/i,T) = yP + a,4 = ^0+5.4 = ^ ln£^(M) (4) 
= 2 + £ In ( tanh - J = 4 |r| + o(r) , r = 1 - t/t c . (5) 



Therefore, according to the relation Eq. ([T]) and Eq. Q with z — 1/3 for two-dimensional 
case this leads to the following result for the surface tension of the liquid- vapor interface of 
the 2D LJ fluid: 

a LJ (T) = *(t(T)) t ^=_ q f (1 - T/T c ) + . . . , (6) 

The result of the surface tension for the 2D LJ fluid is on Fig. [TJ The comparison with 
the result of the available simulational data for LJ fluid is shown on Fig. [2| In the cases 
where no exact solution is available one may use the results of [131 E] • There it was shown 
that the surface tension of the lattice model can be represented as the difference between 
averages of the local physical quantity in coexisting phases: 



° = t ((S) 9 as-(s)u g ) ■ (7) 

The representation ([7]) is based on the markovian property of the distribution function of 
the lattice model with the nearest-neighbor interaction. The distribution function of the 
D-dimensional lattice is represented as the markovian chain for D — 1-dimensional slices. 
Such representation is widely used in the transfer matrix method [15]. Here 5 is the local 
variable which corresponds to the spin distribution in the slice of the gas phase far from the 
interface. In the simplest mean-field approximation of the Bragg- Williams 5 is reduced to: 



0.5 



2D Ising model (Onsager) 



T/T. 



0.2 



0.4 



0.6 



0.8 



1.0 



FIG. 1. The surface tension of the 2D Ising model (Onsager solution) and the surface tension of 
the 2D Lennard- Jones fluid obtained according to Eq. Q. 

where p(si) is the distribution of the z-th spin in the slice. Substitution Eq. ^ into Eq. ([7]) 
and taking into account the symmetry between gas the liquid density: x gas = 1 — X\ iq we 
get: 
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or in terms of the magnetization m — 2x — 1: 



— 1mm I In 1 + m(t) 
1 lJI l-m(t) 



(9) 



(10) 



where Z is the lattice spacing. Further we put l = 1 for simplicity. Sure neglecting the 
correlation leads to the very crude approximation especially in the fluctuational region so 
the direct comparison of Eq. ([7]) with the surface data leads to the big difference [16J. In 
order to correct such weakness but conserving the analytical simplicity for the connection 



between s and x we modify Eq. (10) as following: 
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where £ e ff is the effective width of the interface defined by Eq. (11). It should be of order of 



correlation length. The value a is determined by the correspondence with some reference 
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FIG. 2. The comparison of the surface tension of the 2D Lennard- Jones fluid given by Eq. Q with 
the simulation results of [12] . 



point. For example in case of the 2D Ising model with: 

m(t) = ( 1 1 



sinh 4 ( 2/t ; 



1/8 



(12) 



the value of do is determined by the natural asymptote £/a— )>lat£— >>0 since the thermal 
fluctuations vanish and the interfacial thickness goes to its minimal value. This allows to 



consider Eq. (11) as the definition of the interfacial width. The corresponding result is on 



Fig. [3} Of course such definition relies on the approximation for 5. But we note that near 



the critical point Eq. (11) leads to: 
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(13) 



This has the correct critical asymptote for the surface tension a oc |t|( d-1 ^. With account of 
the global isomorphism relations Eq. ([!]) this transforms to the dependence a oc (nu q — n gas ) 2 
derived in [17] . So the approximation Eq. ^ is sufficiently good to take into account 
fluctuations at least in the leading order. 

Basing on the previous results on the binodals [TD] and knowing that in 3D case the 



fluctuations are less in comparison with 2D case we can use the model Eq. (11) to reproduce 



the data for the surface tension of the LJ fluid from the information of the binodal. Using 
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FIG. 3. Effective interfacial thickness for 2D Ising model from Eq. (11) with ctq = 1/4, 77 = 1/4 on 



the basis of the Onsager's solution Eq. (12). 



the basic relations Eq. ([T]) we get: 



a 



(14) 



Here nu q ^ as are the densities of the coexisting phases along the binodal for which we use 
the Guggenheim binodal expression: 



iq,gas 
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(15) 



4 v ' <" 4 

We also use simple expression for the temperature dependence of the correlation length 
£ e ff —> ^0 = (T c /T — l) _z/ , r] « 0.03 - Fisher's critical exponent and the exponent 1/ is 
taken as the fitting parameter. The result are shown on Fig. [4| We can use this result 



to compare the interfacial thickness determined by Eq. (11) with the simple expression £ 



Such comparison is shown on Fig. [5| The applicability of the expression is restricted by 
that of the Guggenheim's law (valid for 0.55 < T/T c ) and the usage of £o(T) which becomes 
invalid for low temperatures. Indeed the deviation of £ from t; e ff becomes noticeable for 
low temperatures (see Fig. |5(b) ) where the exponential behavior of n on the temperature 
along the binodal must prevail rather the simple polynomial one. To overcome this difficulty 
one can use the results of the simulations for the correlational length of the 3D Ising model. 
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FIG. 4. Surface tension of the 3D LJ fluid, the argon and the fitting according to Eq. (14) and 



Eq. (15) 
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FIG. 5. Temperature dependence of the effective interfacial thickness £ e ff for 3D LJ fluid (a) and 
the ratio ^ e ff to simple expression £o(T) (b). 

As a summary we have demonstrated that the global isomorphism approach based on 
simple transformations Eq. ([T]) can be applied not only for bulk properties of the coexisting 
phases but also for the description of the surface tension. It opens the way to the application 
of the results obtained within the lattice models to the description of the real fluids of the L J 
type. Since the idea behind Eq. ([T]) is heavily based on the restoration of the particle-hole 
symmetry the application to the Tolman length [18] is of special interest. This is due to 
the fact that this finite size correction to the surface tension vanishes for the symmetrical 



7 



models [19J. Therefore it is natural that at least phenomenologically the Tolman length 
can be directly connected with the density diameter [20]. Within the global isomorphism 
approach this means that it is possible to relate the Tolman length with the asymmetry 
parameter a. This will be the subject of the future work. 
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